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THE EFFECT OF INTERNAL PRESSURE ON THE BUCKLING STRESS 
OP TKIN-WALLED CIRCULAR CYLINDERS UNDER TORSION 
By Earold Crate, S. B. Batdorf, and George V/. Baab- 



The results of a s'evxe^ of tests to determine the 
effect of internal pressure on the buckling load of a 
thin cylinder under an applied torque indicated that 
internal pressure raises the shear buckling stress. 
The experimental results were analyzed with the aid of 
previously developed theory and a simple interaction 
formula was derived. 



The curved metal skin of a modern airplane in 
flight is subject to stresses that m.ay cause the skin 
to buckle, and proper design of airplane structures 
requires a knowledge of the stress conditions under 
v/hich buckling will occur. The ability to estim.ate the 
buckling point under combined loading conditions is of 
particular importance . 

In order to determine the effect that norm.al pres- 
sure, or air loads, might have on the critical stress 
for curved sheet, two preliminary tests were miade (ref- 
ences 1 and 2). A pronounced increase in critical 
stress with increase in normal pressure was found, and 
the subsequent interest shown by the aircraft industry 
in this subject indicated the desirability of further 
study. 

No well-established theories for buckling of curx^ed- 
sheet panels either in torsion or under hydrostatic 
pressure are available; however, satisfactory theories 
for buckling of com.plete cylinders under these loadings 
have been advanced (references 3 and 4). In order to 
effect som.e correlation between theory and experiment. 



SUMMARY 
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therefore, an Inves ti^-^ation wa.s made of the influence of 
Internal pressure on the critical stresses of thin 
cylinders in torsion. Experirients were conducted to 
determine the critical shear stresses for four cylinder 
lengths at a number of different internal pressures. 
The' theories and the experimental data were used in con- 
junction to determine an interaction formula for the 
buckling of cylinders of moderate len.^th under the com- 
bined action of torsion and internal pressure. 



SYJffiOLS 

E Young's modulus of elasticity, psi 

G shear modulus of elasticity, psi 

length of cylinder v:ithout rings, inches 

L length of cylinder between rings, inches 

t thickness of cylinder v^all, inches 

d diam.eter of cylinder. Inches 

r radius of cylinder, inches 

p internal pressure of cylinder, considered 

positive v/hen it produces tens'j le stresses 
in cylinder walls, psi 

^ critical pressure in the absence of torsion 

^ cr jrro (negative according to sign convention 

adopted for p) , psi 

T shear stress in cylinder walls due to applied 

torque, psi 

T value of shear stress v/hen cylinder is at the 

point of elastic instability, psi 

) critical shear stress in the absence of 

^ P^'^ interrxal pressure, psi 

Rp pressure ratio 
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Ro shear-stresj^ ratio 



q exponent of Rg in the interaction forrr.ula 

9 rotation of free end of cylinder, radians 

Kg coefficient used in Lvndquist's empirical 

forriula (appendix A) 

y,y^ differences in strain gage readings 

[i Poisson^s ratio 



SPECir^EITS AIID T!:.3T PROCEDUI^vB 



Diagrams of the test-cylinder construction and 
ring systems used are given in figure 1. The cylinder 
v;as made from 0.032-inch 24S-T aluminum-alloy sheet 
closely riveted around tv/o heavy steel rings, one at 
each end. The sheet was joined along an element of the 
cylinder v/ith a butt joint covered with a single strap 
on the outside. Rings made of 24S-T aluminum alloy 
were added to this cylinder and divided it into shorter 
cylinders, the lengths of which were equal to the ring 
spacings . 

Figure 2 is a photograph of the apparatus used to 
test the cylinder. The ends of the cylinder v;ere 
closed by heav^;- flat steel plates in order that air 
could be maintained under pressure inside the cylinder. 
The weights of the steel plate and test ring at the 
free end of the cylinder were neutralized by an upward 
load on one of the torque arms. Rotations of the free 
end of the cylinder relative to the floor were measured 
by a pair of dial gages. 

Observed buckling loads at zero and at low internal 
pressure were determined as the leads at wliich there 
was a sudden snap of the cylinder to the buckled state. 
This snap action v/as in many cases preceded by a slow 
growth of visible wrinkles in the cylinder v/alls as the 
load increased and a greatly increased rate of growth 
of the wrinkles at close proximity to the snap-buckling 
load. With increases in the internal pressure, the 
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5^nap decreased in violence lontil It was no longer 
observable o In these cases the buckling load was esti- 
mated, on the basis of visual observations, to be that 
load at which the rate of growth of wrinkles Y./ith load 
was comparable with the rate of growth of v/r inkles just 
prior to buckling in those cases in which snap buckling 
did oceur. 

The Southv;ell method of determining the critical 
stress of specimens vrith initial eccentricities was 
assumed to be applicable to cylinders subjected to 
torsion and was used as a check on the visual determi- 
nation. Two electrical strain gages were mounted on 
opposite sides of the cylinder wall at the crest of a 
buckle, the location of" which was determined by means of 
a preliminary buckling test. The difference y in 
readings of the strain gages on opposite sides of the 
cylinder wall provided a numerical m.easure of the dis- 
tortion as the wrinkles grew. The Southwell method was 
applied by plotting ^y - ^lyC'^ *" '^l) ^f^^inst y - j-^ 
where y^ and were arbitrarily chosen initial 

values of each quantity (reference 5) . The inverse 

slope of the straight line forr^ied in this r.anner is 

T - where r is the desired critical stress, 

cr 1 cr 

DISCUS3I0i^ OF R^jSTJLTS 



In figure 3 the shear stress calculated from, the 
external load on the cylinder is plotted against the 
corresponding rotation of the free end for various 
values of ring spacing and internal pressure. The 
solid lines give rotations computed by the formula 




in which G v:as assumed to be 3,970,000 psi. 

Prior to buckling, the stiffness of the cylinder in 
torsion (slope of the curves) is essentially unaffected 
by the internal pres'^ure or by the ring spacing. The 
rotations measured were slightly greater than those 
predicted by the formula. This' result may be attributed 
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to the method of neasurenent^ which considered as rota- 
tions many other small effects, such as bending of the 
column supporting the cylinder, distortion of rivets and 
rivet holes, and taking' up slack in the bolt holes in 
the end plates . 

The experimental buckling stress in torsion without 
pressure is com^Rved with the predictions of Donnell 
(reference 3) and Lundquist (reference 6) in the fol- 
lowing table? 





Critical she 


ar stress, 


ksi 


L/d 
















ExperlTiental 


Lundquis t 


Donnell 


0.18 


5.8 


6.0 


7.9 


.36 


3.9 


4.2 


4.9 


.72 


3.0 


3.0 


% % 


1.43 


2.3 


2.3 


2.3 



LiJindquis t ' s formula was derived for stress at failure. 
Except for very short cylinders, however, the shear 
buckling stress and the stress at failure are essen- 
tially the sam.e. Donnell 's theoretical curve (solid 
curve' of fig. 4) and L-ondquis t ^ s empirical form.ula are 
discussed more fully in appendix A. 

In figure 5 the buckling stresses under torsion 
are plotted against the internal pressure for each 
cylinder; this figure indicates tl.iat the buckling 
stresses increase as pressure increases* The results 
were essentially the same whether the visual or the 
Southwell method, of determining buckling loads was 
used, although the Southv/ell method usually gave 
slightly higher loads than the visual. 

The buckled cylinder could be returned to the 
unbuckled state either by decreasing the torque or by 
increasing the pressure. The loads at v/hich the buckles 
disappeared were determined visually for cylinders 1, 
2, and 3a, and are also shown in figure 5. These loads 
were not recorded, however, for the remaining cylinders 
because of a large scatter in the readings. This 
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scatter is attributable to the fact that, as the cylin- 
ders became shorter, the disappearance of the buckles 
became more gradual and the visual selection of a 
"point of disappearance" became difficult* In the case 
of cylinder 2, in which the readings are relatively 
definite, the location of the curve indicating disap- 
pearance of buckles was the same for either an increase 
of pressure or a reduction of torsion. 

The method used to find the interaction formula 
best representing the experimental data is given in 
appendix A. The analysis leads to a formula of the 
tyipe 

R ^ + R =1 
s p 

where the value of the exponent q depends upon the 
ass^am.ptions r.ade for (PcrX- o ^^'^ (^^^)r)-0* "^^^ 
cylinders of mioderate length, which according to 
Donnell's analysis satisfy the inequality 



30* <(hf< 4 

the exponent q is equal to 1-89 to 2 using t^^ 

derived from Donnell^s theoretical curve, the exact 

value depending on the value of iF'/td concerned. 
When T^^ derived from. L^jindquis t » s empirical forriula 

is used, q - 2 • 17 • 

In figure 6, curves representing the exponents 
1.89, 2, and 2.17 are drawTi through the experim.ental 
data of figure 5 replotted by a method explained in 
appendix A. The three curves lie close together and, 
so far as fit of data is concerned, little basis exists 
for a choice among them.. Simplicity and proximity to 
the average value, however, recommend use of q ~ 2. 
The equation then becomes 
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The parameter L /td, which according to theory 
determines buckling behavior, may be varied by changing 
the length, thickness, diameter, or any combination of 
dimensions. For this reason, the restriction of the 
tests to one thickness and one diameter probably does 
not constitute a significant loss in generality. 

The question of the applicability of the form.ula to 
curved panels is discussed briefly in appendix B. 



CONCLUSION 



The critical stress of a cylinder in torsion 
increases as the internal pressure Increases. The fol- 
lowing Interaction formula was found to represent 
approximately the buckling of a cylinder of moderate 
length under the combined effects of torsion and internal 
pressure : 




where R^ is the ratio of critical shear stress with 

internal pressure to critical shear stress without 
Internal pressure and R-^ is the ratio of internal 

pressure to the critical pressure in the absence of 
torsion. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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APPENDIX A 



DERIVATION OF INTERACTION FORMULA 



Introductory D3.s cuss ion 

Any attempt to detemine which of the conventional 
types of interaction formulas best fits the experimental 
data presented in this paper is complicated by the 
question of v/hich values to use for {P^^\ and 

(jcr) ^^^'P^i'-^'^^^^tal values for ('^cr)pziO '^''^^^^ 

obtained^ but the correrponding experimental values of 
/p \ were not investigated because of the destruc- 

^^^cryTrrO 

tive effects on the cylinder involved in such a test. 
A simple experimental result makes possible a derivation 
of an interaction formula without necessitating assump- 
tions concorning stress ratios or the type of formula 
to be used. 



The experimental data appear to indicate that the 
interaction curves for the four cylinders tested are 
identical, differing only in position. (See fig. 5.) 
The effect of a change of length appears to be a shift 
of the curve parallel to the pressure axis. In fig- 
ure 6 this effect is made clearer by plotting all the 
data to a common p-intercept. The curve for ^ = 1,43 

is shifted parallel to the P'-axis a distance equal to 
-/p \ (as p:iven bv equation (5) which follov/s) in 

order to make the p-intercept zeroo Each of the other 
curves is then shifted the proper distance for best 
superposition. The relatively slight scatter shows 
that the experimental curves are nearly superposable • 

If it is assuried that the interaction curves for 
different lengths of cylinder are for practical pur- 
poses identical, a com.putation of the equation for the 
curves from theory is possible. Among the points on 
the curve of figure 6 are four^ indicated by modified 
symbols, representing the special case p = 0 for the 
four values of the length. Because the curve of figure 7 
represents t^^ plotted against P--(Pcr)T=0' ^^^^^ ^^^^ 
points represent (^cr)p^o P"*-^^^^^' against -(Pcr)T=o 
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for various values of the length. The formula for this 
curve can therefore be obtained by elinlnatlng the 
length froTn the equations expressing ('^cr)p-o 

(Pcr)j~o terms of L, and t. The Interaction 

curve for a cylinder of given dimensions can then be 
obtained by shifting the origin to make the p-intercept 
equal to (Per)^-^ ^^^''^ cylinder. The interaction 

form.ula m.ay then be derived by finding sim.plified expres- 
sions for ^T^^^^^^ and (Pcr)T=o' ^li^^^^ti^S 
and transforming to give the final formula. 



Simplified Expressions for ('^cr)p~o 

Two simplified expressions for (j ^^^^^^^ 

derived, the first based on Donnell^s curve (reference 3) 
and the other on Lundquist^? formula (reference 6). 

Borjriell's theoretical result for the shear buckling 
stress of thin-walled cylinders with simply supported 
ends subjected to torsion is given by the solid curve of 
figure 4. For cylinders of m.oderate length, the curve 
is nearly a straight line given by the equation 



•0 • ^3 / ^ X"* 1 • 2G C 



(la) 



while for extremely large values of L^/td a better 
fit is given by the equation of the asymptote 



./L>^-0.5^(^V1.25 



(lb) 



This essentially straifht portion of the curve is 
Included betv/een the limits 



^oi < {if < 4 



(2) 
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At the lower limit of inequality (2), equation (la) ir^ 
in error by about 15 percent and equation (lb), by about 
20 percent. (For the upper li^it see reference 3,) 

Lundquist's forinula is 

v-1.35 



where values of Kc, are given in reference 6 for variou 
values of J./v. Tn order to determine a forr^ula for Kg 
the values of Kg were plotted against L/v on loga- 
rithmic paper. (See fig. 7.) When L/r > 0.5, 




At L/r = 0.55 ( appro.xirpately the lower limit of 
inequality (2) for the cylinder tested), the error in 
this forrrula about 7 percent. Use of this value of 
Kg leads to the forir;ula 

0.46/^W..35 




Simplified Expression for (P^cr^j-o 

A formula based on work by von Mises and developed 
at the David Taylor Model Basin (equation (10) of refer- 
ence 7) for the buckling of a closed cylinder under 
hydrostatic pressure may be written 

/„ ^ = JlSU (4) 

If, as in the inequality (2), 

^4 < (iT 
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the denominator of equation (4) remains within IC percent 
of unlty^ so that, to this degree of approximation, 

The negative sign appears because in the present paper 
internal pressure is considered positive. 



The Interaction Formula 

For the sake of generality and, at the same time, 
of simplicity in the derivation of the interaction 
formula, equations (1) and (5) maybe written 



Ccr)p.o " K2L%(d,t) (6) 



(7) 



where K]_ and K2 arbitrary constant 5:, m and n 

are arbitrary exponents, and f^ Q-nd fo are arbitrary 

functions of diameter and thiclmess, Elim.ination of L 
gives 



7 V -n/n 



1 



This equation is the formula for the curve of fig- 
ure 6, which is a plot of t^^ against p - (Pcr)T-o 
and therefore, v/hen p = 0, a plot of ('^cr) ^S^lnst 
'■(Pcr)j-o' order to find the interaction formula, 

substitute p - (Pcr)T~o pl^^e of ■"(^cr)T~o 

T^^ in place of ('^cr) -q ^'^^ simplify the resulting 
expression. Thus , 
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P - 



Dividing through by ~^Pcr)T-o using equations (6/ 

and (7) yields 

\n/n — _ — ri/n 

Jv-i J- T / i or. 

+ 1 = 



(Per),. 



0 



JVl_/2cr_\ 
K;^L"f^\-2f2/ 



cr 



p=0 



A rearrangement of terns gives 

,rn/n 



Tcr 1 



(^cr) 



P=OJ 



= 1 



0 



If Donnell^s equation (la) is used with equation (5), 



m 



n 0 . 53 



= 1.89 



or if Donnell^s equation (lb) and equation (5) are used. 



n ^ 1 



n Op5 

Use of Lundquist^G equation (S) \vlth equation (5) yields 

= 2.17 



n 



n 0.46 

For purposes of s iriplif ication, let 

T, 



Rs = 



cr 



(^cr) 



p=0 
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and 



m 

q = n 



Then, the final fcrrnula may be written 



Ill- 
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APPENDIX B 



APPLICATIOII OF INTERACTION PORPTOXA TO CURVT.D PANELS 



The theore t?'. cal analysis of the buckling of a 
curved panel is much more difficult than the treat. nent 
of complete cylinders. In the case of buckling under 
torsion alone, the boundary conditions required for 
panels greatly complicate the theoretical analysis. It 
may be presumed, however, that the added restraints 
make the critical stress higher than that of a complete 
cylinder of the same length, thickness, and curvature. 
The relation between the pressure and the direct stresses 
produced in the sheet is also much more complicated than 
in the case of a complete cylinder. The axial stresses 
depend on the area of the end bulkhead and on the areas 
of cross section of sheet and spar elem.ents. It is 
difficult to draw conclusions ref^arding the relation of 
pressure to circumferential stresses. It is therefore 
not to be expected that an interaction formula established 
for com.plete cylinders will necessarily apply to curved 
sheet also. 

Tn order to determine whether the interaction 
formula derived in appendix A may be used to gain a rough 
idea of the strengthening effect of internal pressure in 
the case of curved panels, the formula is applied to the 
test points of the two specimens of reference 2. Figure 
v8 is a sketch of the two specimens tested. In fig- 
ure 9 the experimental shear buckling stresses are 
plotted against the internal pressure. 

Curves representing the interaction formula are 
also plotted in figure 9. In order to apply the formula, 
the experimentally determined values of ('^cr)^rrO ^^^^^ 
used; (Pcr) -g assumed to be the same as for a com- 

plete cylinder of identical length, thickness, and 
curvature, and is found by applying equation ([•.) of 
appendix A. Fair agreement was obtained betv;een the 
curves and the test points and the agreement irlght be 
considerably improved by a more precise method of 



de terrin ning 



(Per) 



T=0 
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Dotted lines represent rings on inside only 



.0335 

J9 




1 



.19 



J 



Cylinders 2 $3a 

.0555 - 



I" 



Cylinders 3bH 
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Figure I. - Test specimens 
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Fig. 2 
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Figure 3.- Shear stress -roTation curves for cylinders 
with various ring spacings and internal pressures. 
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Figure 4.-Crirical shear stress for thin-wallcd cylinders. 
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Figure 5- Effect of internal pressure on critical shear stress for 
thin -walled circular cylinders. 
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Figur« 6.- Datd of figure 5 plotted to adjusted scale and curves derived 
from interaction fornnula for various values of the exponent c| . 
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Figure 7- Value of used in empirical formula for 
critical shear stress of thin-walled cylinders. 
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Section of rib at center line of panel 



Figure 8.-Curved -sheet panels of reference 2. 
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Figure 9.- Interaction formula applied to curved panel speclnnens 
of reference 2. 



